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In a previous work by the authors the one dimensional (doubling) renormalization op- 
• erator was extended to the case of quasi-periodically forced one dimensional maps. The 

theory was used to explain different self-similarity and universality observed numerically in 
the parameter space of the Forced Logistic Maps. The extension proposed was not complete 
in the sense that we assumed a total of four conjectures to be true. In this paper we present 
numerical support for these conjectures. We also discuss the applicability of this theory to 
I the Forced Logistic Map. 
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1 Introduction 

This is the third of a series of papers (together with [22] and [23]) proposing an extension of the 
one dimensional renormalization theory for the case of quasi-periodically forced one dimensional 
maps. These three papers are closely related but each of them has been written to be readable 
independently. See also [20] for a more detailed discussion. In the previous two papers we were 
concerned with the theoretical part of the theory. In this paper we include different numerical 
computations which support the conjectures introduced for the developing of this theory. To 
do that we briefly review the theory developed in the previous two papers, skipping technical 
details and we adding some numerical computations to the discussion. 

The universality and self-renormalization properties in the cascade of period doubling bifurca- 
tions of families of unimodal maps is a well known phenomenon. The paradigmatic example is 
the Logistic map /^(x) = ax(l — x). Given a typical one parametric family of unimodal maps 
one observes numerically that there exists a sequence of parameter values {dnjneN C / 
such that the attracting periodic orbit of the map undergoes a period doubling bifurcation. 
Between one period doubling and the next one there exists also a parameter value s^, for which 
the critical point of is a periodic orbit with period 2". One can observe numerically that 

lim ^'^"^"-1 = lim ^""^"-1 = S= 4.66920. ■■■ (1) 

n^oo dnJ^x - dn - Sn 

This convergence indicates a self-similarity on the parameter space of the family. On the other 
hand, the constant d is universal, in the sense that for any family of unimodal maps with a 
quadratic turning point having a cascade of period doubling bifurcations, one obtains the same 
ratio S. 

To explain these phenomena Collet and Treser ( [26] ) and Feigenbaum ([51 [6] ) proposed simul- 
taneously the renormalization operator. Their explanation was based on the existence of a 
hyperbolic fixed point of the operator with suitable properties. The first proof of the exis- 
tence of this point and its hyperbolicity were obtained with numerical assistance |15l [1]. A 
decade later, Sullivan (see [25]) generalized the operator and provided a theoretical proof of the 
hyperbolicity using complex dynamics techniques. See [17\ f3] for extensive summaries on the 
theory. 

In [21] we presented numerical evidences of self similarity and universality for families of quasi- 
periodically forced Logistic maps. These are maps in the cylinder where the dynamics on the 
periodic variable are given by a rigid rotation and the dynamics on the other variable are given 
by the Logistic Map plus a small perturbation which depends on both variables. This kind of 
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AfRne relationship 



Figure 1: Schematic representation of the bifurcations diagram of the Forced Logistic Map, for 
rotation number equal to oj (left) and 2ijj (right). See the text for more details. 

maps have its origins in studies related to the existence of strange non-chaotic attractors (see 
[U El ISl [9l Uni dS]). Let us describe these numerical evidences with more detail. 



1.1 Numerical evidence of self-similarity and universality for quasi-periodic 
forced maps 

Consider {ga,e} {a,£)&JcM? ^ parametric family of quasi-periodic maps in the cylinder T x M 
of the form 

e = e + uo, 

X = ax{l — x) + eh{6, x), 

with UJ a Diophantine number, a and e parameters and h a periodic function with respect to 9 
which can also depend on a and e. Recall that the Logistic map x = ax{l — x) has a complete 
cascade of period doubling bifurcations. As before, let {dn}neN C / denote the parameter values 
where the attracting periodic orbit undergoes a period doubling bifurcation and {s„}„gN C / 
the values for which the critical point of fg,^ is a periodic orbit with period 2". 

In we computed some bifurcation diagrams in terms of the dynamics of the attracting set. We 
have taken into account different properties of the attracting set, as the Lyapunov exponent and, 
in the case of having a periodic invariant curve, its period and its reducibility. The reducibility 
loss of an invariant curve is not a bifurcation in the classical sense that the attracting set of 
the map changes dramatically, only the spectral properties of the transfer operator associated 
to the continuation of that curve does (see [E]). Despite of this, it can be characterized as 
a bifurcation (see definition 2.3 in [H]). The numerical computations in [11] reveal that the 
parameter values for which the invariant curve doubles its period are contained in regions of the 
parameter space where the invariant curve is reducible. These computations also reveal that 
from every parameter value (a,e) = (s„,0) two curves of reducibility loss (of the 2"-periodic 
invariant curve) are born. This situation is sketched in the left panel of figure [!} 

Assume that these two curves can be locally expressed as (s„ + a'j^{uj)e + 0(e^),e) and (s„ + 
(3'„{uj)e +0{e'^) , e) . In [22] we proved that these curves really exist under suitable hypothesis. We 
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have also given explicit expressions of the slopes a'^iuj) and in terms of the quasi-periodic 

renormalization operator (introduced there). We only focus on a'^{uj), but the discussion for 
is completely analogous. 

The slopes can be used for the numerical detection of universality and self-renormalization 

phenomena. If the bifurcation diagram is self-similar by an affine ratio one should have that 
a'j^{cj) / a'^_-^^{uj) converges to a constant. In [2T] we compute numerically this ratios and we show 
that this is not true due to the fact that when the period is doubled, the rotation number of the 
system also is. What we find is that there exists an affine relationship between the bifurcation 
diagram of the family for rotation number oj and the bifurcation diagram of the same family for 
rotation number 2uj. This is sketched in figure [T} 

Concretely, in [21] we observed numerically the following behavior. 

• First numerical observation: the sequence aj^(w)/a^_^(a;) is not convergent in n. But, 
for uj fix, one obtains the same sequence for any family of quasi-periodic forced map like 

with a quasi-periodic forcing of the type h{6, x) = fi{x) cos{9). 

• Second numerical observation: the sequence a'^{u}) / a'^_i{2uj) associated to maps like 
([2]) is convergent in n when we take the quasi-periodic forcing of the type h{9, x) = 
fi{x) cos{9). The limit depends on u and fi. 

• Third numerical observation: the two previous observations are not true when the 
quasi-periodic forcing is of the type hr^{9,x) = fi{x) cos{9) + r//2(x) cos(20) when r] ^ 0. 
But the sequence a^(ci;)/a'^„^(2c<j) associated to the map ^ with /i = /i^ is ry-close to the 
same maps with h = ho 

In |22j we extended the renormalization operator and we obtained explicit expressions of the 
slopes a'jj^{uj) and /3^(w) in terms of the quasi-periodic renormalization operator. This is reviewed 
in section [2| In [23j we give a theoretical explanation to the numerical observations described 
above in terms of the dynamics of the quasi-periodic renormalization operator. This is reviewed 
in section |3j The novelty in this paper is that we present numerical support to the conjectures 
done in |22l 123] . This numerical support is presented after the statement of each of the con- 
jectures, with the exception of conjecture |Aj which is given in section |4] In Appendix |A] we 
describe the numerical approximation used to discretize the renormalization operator and how 
we use it to compute the spectrum of its derivative. 



2 Existence of reducibility loss bifurcations 

Consider a quasi-periodic forced map like 

F : Tx I T X I 

9 + uj \ (3) 



with / G (T X /, /) . To define the renormalization operator it is only necessary that r > 1 , 



but we restrict this study to the analytic case due to technical reasons. Along section 2.1 it is 



not necessary to require uj Diophantine, but it will be necessary in section 2.2 



The definition of the operator is done in a perturbative way, in the sense that it is only applicable 
to maps f{9, x) = g{x) + h{9, x) with g renormalizable in the one dimensional sense and h small. 
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2.1 Definition of the operator and basic properties 
Preliminary notation 

Let W be an open set in the complex plane containing the interval Is = [—1 — 5,1 + 5] and let 
Mp = {z = x + iy&'C such that \y\ < p}. Then consider B = B{Mp,W) the space of functions 
/ : Bp X W ^ C such that: 

1. / is holomorphic in Bp x W and continuous in the closure of Bp x W. 

2. / is real analytic. 

3. / is 1-periodic in the first variable. 

This space endowed with the supremum norm is Banach. 

Let TZTiiW) denote the space of functions / : W — )• C such that are holomorphic in W, continuous 
in the closure of W, and send real number to real numbers. This space is also Banach with the 
supremum norm. 

Consider the operator 

po: B TZniW) 

f{9,x) ^ [ fi9,x)d9. 
Jo 

Let Bq the natural inclusion of Tl7i{W) into B. Then we have that po as a map from B to Bq is 
a projection ((po)^ = Po)- 



Set up of the one dimensional renormalization operator. 

First we give a concrete definition of the one dimensional renormalization operator before ex- 
tending it to the quasi-periodic case. Actually, this is a minor modification of the one given in 

m- 

Given a small value 6 > 0, let Ais denote the subspace of Tl7i{W) formed by the even functions 
■0 which send the interval Is = [—1 — 5,1 + 5] into itself, and such that ip{0) = 1 and xtp'{x) < 
for X 7^ 0. 

Set a = 1^(1), a' = {1 + 5)a and b' = ilj{a'). We can define ViJZs) as the set of -0 S Ms such 
that a < 0, 1 > 6' > -o', and ii{b') < -a'. 

We define the renormalization operator, TZs : V{TZs) — )• Ms as 

7^5(V')(x) = --0 o ^(ax). (5) 
a 

where a = "0(1). 

Note that, given ijj E V{lZs), one needs to ensure that ijj (aW) C W in order to have TZsii^) well 
defined. With this aim, let us consider the following hypothesis. 
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Figure 2: Boundaries of the sets D (g, |) (sohd hne) and (aD (gi |)) (dashed hne). 



HO) There exists an open set W C C containing and a function $ G ;Bsuch that </> = Pq{^) 
is a fixed point of the renormahzation operator TZ^ and such that the closure of both aW 
and (/)(<I>)(aW) is contained in W (with a := <I>(1)). 

In [TBj, Lanford claims that the hypothesis HO is satisfied by the set 



This set is convenient for him because he works in the set of even holomorphic functions. 

In [15] Lanford introduces a discretization of the (one dimensional) renormahzation operator to 
give a computer assisted proof of the contractivity of the operator. In the present paper we use 
the same techniques to discretize the quasi-periodic renormahzation operator, although we do it 
without the use of rigorous interval arithmetics. More details on this are given in the Appendix 
[A} We can use this discretization to check the hypothesis HO for a suitable set W. 

Our study is not restricted to the case of even functions, therefore the set ([6]) used by Lanford 
is not valid in our case. Using the method described in the Appendix [A] we recomputed the 
fixed point of the (one dimensional) renormahzation operator TZ. The fixed point has been 
computed by means of a Newton method with our discretization and then we have checked that 
the Taylor expansion around zero coincides with the one given in [15] . 

With this approximation of (/> we checked (numerically) that D (g, |) the disc of the complex 
plane centered at ^ with radius | satisfies the conditions required to the set W in hypothesis HO. 
In other words, we checked that (j) (aD (g, |)) is contained inside D (g, |) (recall that a = (f){l) 
and aIi{zQ, p) = {z G C | az E ©(zq, p)} ). 

Denote by dn{zQ^p) the boundary of the disk D(zo,/9). In figure [2] we have plotted the sets 
9D(|,|) and (j) [adM which give a visual evidence of the inclusion. Recall that 6 is 

analytic, then to check that the set D (g, |) is mapped inside the set delimited by (j) (a5D (5,5)) 




(6) 
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Definition of the renormafization operator for quasi-periodically forced maps 

Consider the space X dB defined as: 

X = {f C'-iT X h,h)\pdj) Ms}. 

Consider also the decomposition X = Xq® Xq given by the projection pQ. In other words, we 
have Xq = {f G X \ po{f) = /} and Xq = {/ G | po{f) = 0}. Note that from the definition of 
X it follows that Xq is an isomorphic copy of Ms- 

Given a function g & X, we define the quasi-periodic renormalization of g as 

[%{g)]{0,x) := ^g{e + LJ,g{9,ax)), (7) 

where a= g{0, l)d6. 
Jo 

Then we have that there exists a set T>{T), open in (po o TL) ""^ {Ms), wlierc the operator is well 
defined. Moreover this set contains T>q{T), the inclusion of V{JZ) in B. By definition we have 
that %j restricted to I'o(T) is isomorphically conjugate to 7?., therefore the fixed points of Tl 
extend to fixed points of 7^. Assume that HO holds and let $ be the fixed point given by this 
hypothesis. Then we have that there exists U C T^{T) H B, an open neighborhood of <1>, such 
that Tu, '■ U ^ B \s well defined. Moreover we have that Toj is Frechet differentiable for any 



Fourier expansion of DTuj{^). 

Let ^' be a function in a neighborhood of $ (given in hypothesis HO) where 7^ is differentiable. 
Additionally, assume that ^ G T^oiToj). 

Given a function f E B we can consider its complex Fourier expansion in the periodic variable 

f{d,z) = J2ckiz)e^''''\ (8) 

with ^ 

Jo 

Then we have that D%j "diagonalizes" with respect to the complex Fourier expansion, in the 
sense that we have 

[DrU^)f] {9, z) = DUsicoKz) + Yl {[Liickm + [L2(c,)](z)e2-'=-) e^^'^^', (9) 

fcez\{o} 

where 

Li : nn{w) nniw) 

g{z) -ip' o^p{az)g{az), 

and 

L2 : nn{w) nn{w) 

9{z) 1-^ -^goi){az), 



7 



with ip = po(^) and o = ip{l). 

An immediate consequence of this diagonahzation is the following. Consider 

Bk := {/ G B\ f{e,x) = u{x) cos{2TTke) + v{x) sm{2Trke), for some u,v £ nn{W)}, (10) 
then we have that the spaces Bk are invariant by DT{^) for any A; > 0. 
Moreover DTu){^) restricted to Bk is conjugate to Ckuj, where Ca is the defined as 

: nn{w)®nn{w) nn{w)®nn{w) 



u\ I— r "^^(^^ I + ( '^^^i'^'^^) — sin(27ra) \ / L2{u) 

V J y Li{v) J y sin(27ra) cos(27ra) J y L2{v) 



(11) 



Then we have that the understanding of the derivative of the renormalization operator in B is 
equivalent to the study of the operator C^j for a any w G T. 



Properties of £^ 

Given a value 7 G T, consider the rotation defined as 



R^ : n-H{w)enn{w) Tzn(w)enn{w) 



u \ j cos(27r7) — sin(27r7) \ f u 

V J y sin(27r7) cos(27r7) J \ v 



(12) 



then we have that C^j and i?^ commute for any w, 7 G T. 

This has some consequences on the spectrum of Ci^. Concretely, we have that any eigenvalue 
of Cuj (different from zero) is either real with geometric multiplicity even, or a pair of complex 
conjugate eigenvalues. On the other hand C^^ depends analytically on u, which (using theo- 
rems III-6.17 and VII-1.7 of [13j) imply that (as long as the eigenvalues of C^j are simple) the 
eigenvalues and their associated eigenspaces depend analytically on the parameter uj. 

Finally, doing some minor changes on the domain of definition, we can prove the compactness of 
the operator £(^. Recall that the compactness of an operator implies that its spectrum is either 
finite or countable with on its closure (see for instance theorem III-6.26 of jl3j). 

In figure [3] we have a numerical approximation of the spectrum of the operator C^o depending 
on UJ. We can observe that the properties described above are satisfied. The details on the 
numerical computations involved to approximate the spectrum are described in Appendix [A| 
Several numerical tests on the reliability of the results are also included there. 



2.2 Reducibility loss and quasi-periodic renormalization 

Given a map F like ([s]) with f £ B and a; G T, we denote by /" : T x R — )• M the x-projection 
of F"'{x,6). Equivalently /" can be defined through the recurrence 

r{e,x) = f{e + {n-i)u,r~\e,x)), f{9,x)^x. (is) 



8 




Figure 3: Numerical approximation of the spectrum of for w G T. Top: projection in the 
complex plane of the spectrum when u varies in T. Bottom left: evolution of the real part with 
respect to w. Bottom right: evolution of the imaginary part with respect to oj. 



From this point on, whenever u is used, it is assumed to be Diophantine. Denote hy Vt = Vt-y^r 
the set of Diophantine numbers, this is the set of a; G T such that there exist 7 > and r > 1 
such that 

\qoJ-p\>^, for ah G Z X (Z\{0}). 

Additionally, we will need to assume that the following conjecture is true. 

Conjecture A. The operator Tuj (for any oo ^ Vt) is an injective function when restricted to 
the domain Br\V{T). Moreover, there exists U an open set of V{T) containing W^i^^^TZ) U 
Vl^''(<I', 7^1^ if/iere the operator Tuj is differentiahle. 

In |22j we discuss the difficulties for proving this conjecture. A priory there is no way to check 
numerically this kind of conjecture. A posteriori we have that the results obtained assuming 
this conjecture are coherent with the numerical computations (see section |4]). 

^Here W°{'^,TZ) and W"("l>,7?.) are considered as the inclusion in B of the stable and the unstable manifolds 
of the fixed point $ (given by HO) by the map TZ in the topology of Bq (the inclusion of one parametric maps in 
B). 
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Consequences for a two parametric family of maps 

Consider a two parametric families of maps {c(q, e)}(^a,,s)eA contained in B, with A = [a, b] x [0, c] 
and a, b and c are real numbers (with a < b and < c). We assume that the dependency on the 
parameters is analytic. 

Consider the following hypothesis on the family of maps. 

HI) The family {c{a, £)}(a,e)&A uncouples for e = 0, in the sense that the family {c{a, 0)}a^[a,b] 
does not depend on 6 and it has a full cascade of period doubling bifurcations. We assume 
that the family {c(a, 0)}a^^a,b] crosses transversally the stable manifold of the fixed point 
of the renormalization operator, and each of the manifolds S„ for any ti ^ 1 , where X]^ is 
the inclusion in B of the set of one dimensional unimodal maps with a super-attracting 2" 
periodic orbit. 

In other words, we assume that the family c(a, e) can be written as, 

c{a, e) = co{a) + eci(a, e), 

with {co{a)}a^[a,b] ^0 having a full cascade of period doubling bifurcations. 

Given a family {c{a,£)}(^a,e)€A satisfying the hypothesis HI, let a„ be the parameter value for 
which the uncoupled family {c{a,0)}a^[a,b] intersects the manifold Note that the critical 
point of the map c{an,0) is a 2"-periodic orbit. Our main achievement in [22] is to prove that 
from every parameter value {an,0) there are born two curves in the parameter space, each of 
them corresponding to a reducibility loss bifurcation. If we want to give a more precise statement 
of the result we need now to introduce some technical definitions. 

Let TZ7i(Ep,W) denote the space of periodic real analytic maps from Bp to W and continuous 
in the closure of Bp. Consider a map fo £ B and u € Q,, such that / has a periodic invariant 
curve xq of rotation number uj with a Lyapunov exponent less equal than certain —Kq < 0. 
Using lemma 3.6 in j22j we have that there exist a neighborhood V G B oi fo and a map 
X G TZT-L{Mp,W) such that x(/) is a periodic invariant curve of / for any f G V. Then we can 
define the map Gi as 

Gi : QxV 7^7^(Bp,C) 

{u,g) ^ D,g{e + cj,9{e,[x{u,g)]{9)))D,g{e,[x{u,g)]{e)). ^^^^ 

On the other hand, we can consider the counterpart of the map Gi in the uncoupled case. Given 
a map fo G Bq, consider U C Bq a neighborhood of /o in the Bq topology. Assume that /o has a 
attracting 2-periodic orbit xq G I. Let x = x{f) G W be the continuation of this periodic orbit 
for any f £ U. We have that x depends analytically on the map, therefore it induces a map 
X : {/ — 7- W. Then if we take U small enough we have an analytic map x : U ^ W such that 
x[f] is a periodic orbit of period 2. Now we can consider the map 

Gi: UcBo ^ C 

(15 1 

/ ^ Dj{f{x[f]))Dj{x[f]). 

Note that Gi corresponds to Gi restricted to the space Bq (but then Gi (/) has to be seen as an 
element of nn{Mp,W)). 
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Consider the sequences 

Uk = 2ujk-i, for /c = 1, n — 1. 

^(n) ^ 7^(/i!\), for A: = l,...,n-1. 

4") = Dn(ft\)u^^\, fork = l,...,n-l. 



(16) 



4"^ = ^7L,_, (/S) forfc = l,...,n-l. 



with 



= c(a„,0), m[,"^ = 9„c(an,0), w^"^ = 5£c(a„, 0). (17) 

Note that /^"^ = {c(a,0)} n S„, then /(J"^ tends to W'{n,^) when n grow. Therefore, the 
sequence {ft^}o<k<n attains to W'{n,^) U VF"(7^,$) when n grows and consequently there 
exist no such that {/^"^}o<fe<n C U , where U is the neighborhood given in conjecture [a[ If the 
conjecture is true, then the operator is differentiable in the orbit {/^"^}o<fe<n C U. 

Consider the foUowing hypothesis. 

H2) The family {c{a,e)}(^a,e)&A is such that 

DGi Ln-i, &) D%^_, (/tU • • • D%, 9,c(a„, 0), 



has a unique non-degenerate minimum (respectively maximum) as a function from T to 
M, for any n > tiq. 

Consider a family of maps {c{a,£)}(^a,e)&A such that the hypotheses HI and H2 are satisfied 
and ujQ £ If the conjecture A is true, then theorem 3.8 in [22J asserts that there exists no 
such that, for any n > uq, there exist two bifurcation curves around the parameter value (a„, 0), 
such that they correspond to a reducibility-loss bifurcation of the 2^-periodic invariant curve. 
Moreover, these curves are locally expressed as {an + a'^{i^)£ + o{£),£) and {a~+l3'^{uj)e + o{e),£) 
with 

= ^ ^- ^ 1 ^, (18) 



and 



PM = ^ ^ (19) 

DG, 



where Gi and Gi are given by equations (14) and (15), and m and M are the minimum and the 
maximum as operators, that is 

m: nn{Mf„c) K 

■ m\ (20) 
q I—)- mmg c^i. 

and 

M: nn{Mf„c) K 

m\ (21) 
q I—)- max(7(y). 
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Let us focus again on hypothesis H2, which is not intuitive. We can introduce a stronger 
condition which much more easy to check. Moreover this conditions is automatically satisfied 
by maps like the Forced Logistic Map. Consider a family of maps {c{a,e)}^a,e)eA before, 
satisfying hypothesis HI. 



H2') The family {c{a,e)}ta,e)eA is such that the quasi-periodic perturbation di;c{a,0) belongs 



to the set Bi (see equation (10)) for any value of a (with (a, 0) G A). 



Then we have that H2' implies H2 (see proposition 3.10 in |22j). 

3 Asymptotic behavior of reducibility loss bifurcations 

Consider a two parametric family of maps {c(a, e)}(a,e)eA contained in B, with A = [a, b] x [0, d] 
and o, b and d are real numbers (with a < b and < d). We assume that the dependency on 
the parameters is analytic and the family is such that the hypotheses HI and H2 introduced 



in section 2.2 are satisfied. Consider also the reducibility loss bifurcation curves associated to 
the 2^-periodic orbit given by ( |18| ). Since the value a'^ioo) depends also on the family of maps c 
considered, we will denote it by a'^iuj, c) from now on. We omit the case concerning /3^(u;) (see 



(19)) because it is completely analogous to the one considered here. 

3.1 Rotational symmetry reduction 

Given 7 G T, consider the following auxiliary function 

B 

(22) 



t^: B ^ B 



v{9,z) I-?- v{9 + ^,z). 



Let Bi be the subspace of B defined by (10) for A; = 1. The space Bi is indeed the image of the 
projection tti : B ^ B defined as 

[ni{v)] {6, x) = v{e, x) cos{2TTx)de^ cos(27r6l) + v{9, x) sm{27rx)d9^ sin(27r6l). (23) 

Given xq G W n M and 9o £T we can also consider the sets 

B[ = B[{9o,xo) = {feB, \f{9o,xo) = 0,def{9o,xo) > 0}, 

and 

B' = B'{9o,xo) = {f € B\7T,{f) B[}. 

Note that B'i{9o,xo) depends on the election of {9q,xo), but for any fixed xq and the set 
B'i{9o,xo) is an open subset of a codimension one linear subspace of Bi. Note also that any v G 
then V ^ due to the condition dgv{9Q,xo) > 0. Moreover for any v £ Bi\ {0} there exists a 
unique 70 G T such that t^o{v) G B'i{9o,xo) (see proposition 3.1 in [23]). 

Consider a two parametric family of maps {c(a, e)}i^a,£)eA contained in B satisfying the hypothe- 



ses HI and H2 as in section 2.2 Consider also the reducibility loss bifurcation curves associated 
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to the 2^-periodic orbit with slopes given by (18) and (19). Then we have that the formulas 



(18) and (19) can be expressed in term of vectors in B'i{9o,xo). Let us see this with more detail. 
Consider the sequences {w^}, {/^"^} and {w^"^} given by (16) and (17). Consider now the 



sequence 



and 



4"^ = ho {dec{an,0)) 



(24) 
(25) 



where tI'V^,"^) and 70 are chosen such that v]^' belongs to B'{6o,xq) for A; = 0, 1, ...,n. 



If the projection of DTu]^_^ (^fk^ij ^k^-i in tSi is non zero, then 7 is uniquely determined 

(n) 

and the vectors vj^ are well defined. Moreover, if we assume that ujq G 0, and that the conjecture 



A is true then we have that a'j^{uj) can be written as 



(26) 



where Gi, Gi and m are given by equations (14), (15) and (20). For more details see theorem 
3.2 in p3]. 



3.2 Reduction to the dynamics of the renormalization operator 

The goal of this section is to reduce the problem of describing the asymptotic behavior of 
a'j^{u)o, ci) / a'^_i{uJo, ci) to the dynamics of the quasi-periodic renormalization operator. 

Definition 3.1. Given two sequences {ri}i,zz^ and {sjjjg^^ in a Banach space, we say that 
they are asymptotically equivalent if there exists < p < 1 and ko such that 

W^i — Sill < kop^ Vi G Z+. 

We will commit an abuse of notation and denote this equivalence relation by ~ instead of 

Given a family {c{a,e)}(^a,e)&A satisfying hypotheses HI and H2 and a fixed Diophantine rota- 
tion number uq, let a* denote the parameter value such that the family {c(a, 0)}(q, o)eA intersects 
with W^**($,7^) and denote the intersection of W'"'{^,TZ) with the manifold Sj. Consider then 



ojk = 2cL!fc_i, for A; = 1, n — 1, 

f DTZ ($) Uk-i, for A; = 1, ... , [n/2] - 1, 
\ DTZ Uk-i, for k= [n/2],..., n-1. 



(27) 



Vk 



with 



t^^^^_^) {D%,_, ($) Vk-i) , for fc = 1, [n/2] - 1, 

h{vk-i) (^"T^ft-i (fn-k) Vk-i) , ior k= [n/2], ...,n-l. 

uo = dac{a*,0), VQ = (9ec(a*,0)) , 
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n 



4 
5 
6 
7 
8 
9 

10 
11 
12 
13 



n/2]-l\ 



5.426626e+01 
5.426626e+01 
2.361437e+02 
2.361437e+02 
1.158124e+03 
1.158124e+03 
6.859354e+03 
6.859354e+03 
5.625187e+04 
5.625187e+04 



m 



DGa 



1 I ^[n/2]-li 



•"[.1/2]-! 
lk[n/2]-lll 



1.666220e+01 
1.666220e+01 
1.394208e+02 
1.394208e+02 
5.923824e+02 
5.923824e+02 
2.901964e+03 
2.901964e+03 
1.717016e+04 
1.717016e+04 



-{") 



\'"[n/2]-l I 



9.798573e-02 
1.019194e-01 
4.224092e-03 
4.246578e-03 
1.410327e-03 
1.381662e-03 
4.375003e-04 
4.174722e-04 
1.071136e-04 
1.024221e-04 



Table 1: Different values related with the conjecture B for the Forced Logistic Map (37) and 

, , _ 

^0 - ^■ 



and '^{vs-i) and 70 are chosen such that Us""^ belongs to B'i{9o,Xo) for any s = 1, 
Consider the following conjecture for the forthcoming discussion. 

Conjecture B. For any family of maps {c{a,e)}(^a,e)&A satisfying HI and H2, assume that 



-in) 



I ~(") II 117) 1 I 

In— ill 



with v^li and Vn-i given by (24) and (21). Also assume that there exists a constant C > such 
that 

||f„_i|| > C for any n > 0. 
Finally assume that there exists a constant Co > such that 



m DGi Un^i,fl, 



Vn-l 
\Vn-l\ 



>Co, 



for any n > anduQ Diophantine, where m is given by (20), Gi by (I4) and {fi} = W'^{TZ, $)n 
Si. 



(n) 

In the first part of the conjecture we assume that the asymptotic behavior of the vectors w„_i 
is determined by the linearization of the renormalization operator in the fixed point. We have 
that the iterates of /g"''' correspond to a passage near a saddle point. The initial point /q"'* is 

always in {c(a, 0)}(q o)eA; the final point f^i is always in Ei for any n, and the orbit of the 
points spends more and more iterates in a neighborhood of <1> when n is increased. Therefore 
it is reasonable to expect that the asymptotic behavior is governed by the linearization of the 
operator on the fixed point. In the second part of the conjecture we assume that the modulus 
of the vector does not decrease to zero. 

Table [1] supports numerically conjecture B for the Forced Logistic Map (37) with lo = Note 
that instead of computing the values ||vn-i 
which appear in conjecture B 



we have computed ||t'[„/2]-il|) ""^ [DGi (i^[„/2 



and 



(n) 

n-1 fn-1 



^ ^[n/2]-l 

-1'^' \hnm-i\\JJ 



14 



-(") 



and II — II '"'^^^ ^11 ■ We have done this basicahy to avoid computing the maps /* = 

W^i^^ ^TV) n Sj, which would require computing 1^"(<I>,7^) the unstable manifold. The point /* 
accumulate to the fixed point $ when i — )• oo, then one should expect the same behavior of the 
sequences which appear in conjecture B and the sequence computed in table [T] For more details 
on how this values are computed see section [4j 

Finally we will need the following extension of hypothesis H2 



H3) Consider a two parametric family of maps {c(a, e)}(a,e)eA (with A C M^) satisfying HI and 
H2 and a fixed Diophantine rotation number uj^. Consider also a;„ and Vn given by (27) 

and the point {/i} = VF"(7^, n Si. We assume that DG\{ijJn-\-, ^\)v'^}_y has a unique 
non-degenerate minimum for any c jq £ and n > 0. Assume also that the projection of 

' „-,(") 



An) 
Jk-1 



\ in Bi given by 
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IS non zero. 



Let {c{a,e)}(^a,e)&A (with A C M^) be a two parametric family of q.p. forced maps satisfying 
HI, H2 and H3 and let ojq be a Diophantine number. Consider the reducibility-loss directions 
a'^{c,ijJo) and the sequences Un and Vn given by (27). In theorem 3.6 of [23] we proved that, if 



conjectures A and B are true, then 



a'„_i(c,wo) 



m I DGi { LJn-i, n 



Vn-l 
\Vn-l\ 



m DGi ujn-2,fl 



Vn-2 
\Vn-2\ 



Vn-2 



\Vn-2\ 



(28) 



where m is given by d by ([14|), {f^} = W^iTZ, <i>) n Si is the intersection of the unstable 
manifold of TZ at the fixed point ^ with the manifold Si and S is the universal Feigenbaum 
constant. 

This reduces the asymptotic study of the ratios to the sequence of vectors , which 

are determined by the dynamics of the q.p. renormalization operator. 



3.3 Theoretical explanation to the first numerical observation 

Consider a two parametric family of maps {c(a, s)}[a,e)&A contained in B satisfying the hypothe- 
ses HI, H2 and H3. Let ujq be a Diophantine rotation number for the family. 



The values f"*-*^"'^^ depend only on the sequences Un and Vn given by (27), with vq = dsc{a* , 0) 



-i(^o,c) 

and a* the parameter value for which the family intersects W^{TZ,^). T. 
Vn is described by the dynamics of the following operator. 



le behavior of vectors 



L: T xB' ^ T xB' 

t^^,)iDUm \ (29) 



where 7 is chosen such that t^(^) {DTuj{^)v) belongs to B' . 



In order to study numerically the map L let us recall that the spaces B^ given by (10) are 



invariant by DTw{^), moreover the restriction of D7lo{^) to this space is equivalent to the map 
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C-kuj given by (11) with the space Bi identified with TZTiiW) © TZTi{W). For more details see 



proposition 2.16 in |22| . 

Let us consider the following maps: 

Li: T xBi ^ T xBi 



and 

Li: TxB[ ^ T X B[ 



10 ^'"^ ^ 
(w, v) I— 7- I zo;, 



(31) 



with 



(32) 

where ^{v) is chosen such that t^(^) {Ci^{v)) G fi'^^. 

We have that L restricted to B'l is equivalent to the map Li. When L is restricted to Bk with 
k ^ 1, then it is equivalent to Li. Actually we have that the map Li is equi valent to the map 
Li after applying the rotational symmetry reduction described in section 



3.1 



We can use the discretization of £^ described in Appendix [a| to study numerically the maps Li 
and Li. Let us focus first on the case concerning Li. Given v G Bi = Tn-l{W)QTl'H{W), consider 
the coordinates of v = {x, y) given by this splitting. Following the discretization described in 
Appendix [A] each function x G TZT-LiW) can be approximated by the vector {xo,xi,X2, ...,xn) G 
M^+i where Xi is the i-th coefficient of the Taylor expansion for x around 0. This also holds for 
y, the second component of v. With this procedure each element v in Bi can be approximated 
by a vector {xo,xi, ...,XN,yo, ■■■,yN) in M2(^+^^ 

Following the same argument we can approximate Li : T x Bi — t- T x ;Bi by a map : 
T X M2(^+i) — ^ T X M2(^+^). We can use the discretized map to study the dynamics of 
Li. Given an initial point vq = {x^,y^) ^ 0, we have iterated the point by the map a certain 
transient A^i and then we have plotted the following N2 iterates. Figures |4] and [5] show different 
projections of the attracting set. The values taken for this discretization are N = 30, A^i = 2000 
and N2 = 80000. We have displayed the coordinates corresponding to the first even Taylor 
coefficients of the functions x and y. The odd Taylor coefficients obtained were all equal to zero. 
This last observation suggest that the attractor is contained in the set of even functions (note 
that the subspace of Bi consisting of all the even functions is invariant by C^)- 

The same computations have been done for bigger values of and the results are the same. 
This indicates that the set obtained is stable with respect to the order of discretization, therefore 
it is reasonable to expect that it is close to the true attracting set of the original system. 

Let us remark that we have not made explicit the initial values of wq and vq taken for the 
computations. Indeed, the results seem to be independent of these values. We have repeated 
this computation taking as initial value of vq all the elements of the canonical base of the 
discretized space M^(-'^+^) and we have always obtained the same results. We have also repeated 
the computations for several values of ojq obtaining always the same results. 

Given a solid torus T x Dp, with Dp the disk of radius p in C, we have that the map f{9,x + iy) = 
{cos{9){x + Kq) , sin{6){x + Kq) , y) embeds this torus in M^, for any Kq > p. This embedding can 
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Figure 4: Several planar projections of the attractor of the map (|30j). Form left to right, and top 
to bottom we have the projections in the coordinates (w,xo), {uj,yo), {xo,yo), (w,X2), {u;,y2), 
(3:2,2/2), (w,X4), (w,2/4) and {xi,yi). 



be used for a better visualization of the spatial projections of the attracting set of the map (31 ). 
In the right hand side of figure [5] we have plotted the image by the embedding of the points on 
the left hand side. The concrete values of Kq taken are (from above to below) 2, 3/4 and 3/20. 

The numerical approximation of the attractor displayed in figures |4] and [5] reveal the rotational 
symmetry of the attractor. This is the rotational symmetry described in section [3Tj 



We can use the discretization described in Appendix|X]to approximate numerically the dynamics 
of Li, as we have just done in the case of Li. For the numerical simulation of the operator, we 
have taken = and xq = 0. In this case the set B[{0,0) is identified in M^CiV+i) ^ith the half 
hyper plane 

{(x, y) £ m2(^+i) I xo = and yo > 0}, 
where xq and yo are respectively the first components of x and y. 

In figures |6] and [7] there are displayed different projections of the attracting set obtained iterating 
the map Li. As before we have also considered the map f{9,x + iy) = {cos{9){x + KQ),sm{6){x + 
Ko),y) which embeds the solid torus T x Dp in for a better visualization of the set. This 
time the values of Kq have been taken equal to 1/2 (above) and 3/25 (below). 

Note that the different projections of the attracting set displayed in figure [6] keep a big resem- 
blance with the plots of the dyadic solenoid displayed in figure 5 of [18] . Indeed we believe that 
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Figure 5: Several spatial projections of the attractor of the map (30). In the left hand side of 
the picture we have a plot (from top to bottom) of the projections in the coordinates (to, xo, yo), 
{uj,X2,y2) and (tj, 2:4, 7/4). In the right hand side there are displayed the image of the left side 
projections taking a map that embeds the solid torus in (see the text for more details). 

the attractor is the inclusion of a dyadic solenoid in B[. For more details on the definition and 
the dynamics of the solenoid map see [21 [HI [181 [2l] . To explain this fact, let us introduce a new 
conjecture, this time on the operator Li. 

Conjecture C. There exist an open set V C B[ (independent of u) such that the second 



component of the map Li given by (31) is contractive (with the supremum norm) in the unit 
sphere and it maps the set V into itself for any a; € T. Additionally we will assume that the 
contraction is uniform for any a; € T, in the sense that there exists a constant < p < I such 
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Figure 6: Several planar projections of the section of attractor of the map ( |31[ ). Form left to 
right, and top to bottom we have the projections in the coordinates {u},yo), {uj,X2), {uj, 1/2), 
{X2,y2), {(^,X4), {uj,y4) and (x4,y4). 

that the Lipschitz constant associated to the second component of the map L\ is upper bounded 
by p for any lo GT. 



A good reason to think that conjecture C is true resides in the spectrum of the operator C'^. In 
figure [8] we have a numerical approximation of this spectrum with respect to w as a parameter. 
For this computation we have followed the same procedure that we used for the computation of 
the spectrum of C^), for details see Appendix [A} 

Looking at figure [s] we can observe that there exists a dominant eigenvalue (which is plotted 
in a dashed line) that does not cross the rest, which varies "nicely" with respect to cj. Then 
for each value of w, the normalization of C'^ is a contraction in the sphere, with the eigenvector 
associated to the dominant eigenvalue as a fixed point. This means that conjecture [C] is true 
"point-wise", but this is not enough because the domain of contr activity might depend on u. 



Let us justify now why conjecture C would explain the numerical results obtained for the 



attractor of Li. Consider the set T x y c T x B[{xo,6o). If the conjecture C is true, then 
we would have that the set v invariant by the map Li. More concretely we would have that the 
set would be expanded to the double of its length on the T direction and contracted in the B'^ 
direction. Assume that this transformation is done in such a way that Li maps the set T x V 
inside itself but without self intersections. When we consider the intersection of Li(T x V) with 
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Figure 7: Several spatial projections of the intersection of the attractor of the map ( |31| ). The 
figures in the left correspond to the projection to the coordinates (a;, X2, 2/2) (top) and (oj, x^, y^) 
(bottom). In the right hand side there are displayed the image of the left side projections taking 
a map that embeds the solid torus in (see the text for more details). 



a set of the type {70} x V (for some 70 S T) the section is conformed by two different sets 
without self intersections. The subsequent images by Li we would have (for each leaf {70} x V) 
the double of components than in the previous step, each of them strictly contracted in the B'^ 
component and contained in the previous set. Note that the described process is completely 
analogous to the geometric construction of a dyadic solenoid, but this time contained inside the 



Banach space B'^ instead of the solid torus. Therefore conjecture C would give an explanation 
to the numerical approximation of the attracting set of Li obtained before. 

Consider {ci(a,e)} and {c2(/3,e)} two different families of two parametric maps satisfying hy- 
potheses HI, H2' and H3. Since the family of maps satisfy hypothesis H2', we have that 
9eCj(a*,0) belongs to Bi for i = 1,2. Therefore the dynamics of L (29) coincide with the 



dynamics of Li (31). Theorem 3.10 in |23j asserts the following. 



Consider {ci(a,e)} and {c2(/3,e)} two different families of two parametric maps satisfying hy- 
potheses HI and H2'. Let a* and /?* be the parameter values where each family ci(a,0) and 
C2(/3,0) intersects Ty*(7^, $), the stable manifold of the fixed point of the renormalization op- 
erator. Let Rot(y) = {v ^ Bi \ t^{v) £ V C B'l for some 7 G T} where V is the set given by 



conjecture C for some 7 G T. If 9eCj(a*,0) belongs to Rot(F) for i = 1,2 and the conjectures 
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Figure 8: Numerical approximation of the spectrum of CJ^ with respect to the parameter a;. 
From left to right we have the real part, the imaginary part and the modulus of the first eight 
eigenvalues of L'^ with respect to w. 



A, B and C are true, then (for any € fi) we have that 



a'„_i(wo,ci) a;_i(a;o,C2)^ 



(33) 



where Q;^(a;o,Cj) are the reducibility loss directions associated to each family Cj for the rotation 
number of the system equal to wq . 



3.4 Theoretical explanation to the second numerical observation 

As before consider a two parametric family of maps {c(a, e)}(Q,^e)g^ satisfying hypothesis HI, 
H2' and H3' and let a^(a;,c) denote the slope of one of the curves of the reducibility loss 
bifurcation associated to the 2" periodic invariant curve of the family. In theorem 3.13 of 
23] we prove that the second numerical observation done in section |1.1| can be explained as 



a consequence of the universal behavior (33). One of the hypothesis to prove this requires 
a^(a;, c)/a^(2a;, c) to be a bounded sequence. 

The boundedness of this sequence can be obtained on its turn as a consequence of the following 
conjecture on the operator 



Conjecture D. Consider . 
(11) and ojo G il. Given wo.i o-^'d fo,2 
sequences 



Cca : nniWp) e nn{Wp) nn{Wp) e nn{Wp) the map given by 
'0,1 and fo,2 two vectors in TZTiiWp) © TZ'H{Wp) \ {0}, consider the 

for k = 1, n — 1. 
_i,i) fork = l,...,n-l. (34) 





= 2a;fe_i 


Vk,l 




Vk,2 





Then, there exist constants Ci and C2 such that 

ll'yo,2|| . ||Wn,2|| . ^ ||f0,2| 

^iTi n — 11 n — ^271 T 

lho,l|| |Pn,l|l IFO,iI 

/or any n > 0. 
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Figure 9: We have 



||fO,l|| ||'''n,2|' 
I|f0,2|| ||fn,l|i 



with respect to n, from n = to n = 20 in the left and from 



n 



to n = 2 X 10 in the right. 



This conjecture can be interpreted as a uniform growth condition on C^^. To support this 
conjecture we have computed numerically the iterates J 34 ) , in order to estimate the growth of 

Iko.ill " 



we have plotted the ratios 



\\'"n,2\\ 



\\vn,i\\ with respect to the growth of ||i'n,2||- In figure _ ^ l|t>0 2|| ii^ i| 

with respect to n for the sequence (34) with uo,i = 1 ^-nd vq^2 = 1- For other initial vectors we 
obtain the same behavior. This suggests that conjecture D is true. 



Given {c(a,e)} a two parametric maps satisfying hypotheses HI, H2' and H3 let a* be the 
parameter values for which the family {c(a, 0)}(q o)eA intersects W^iJZ^^^). Let Rot(y) = 



{v ^ Bi \ t^{v) <Z B'l for some 7 G T} where V is the set given by conjecture C for some 



7 G T. Assume that dec{a* , 0) G Rot(y). In corollary 3.14 of [23] we prove that, for any loq G VL^ 
then 

lim (35) 

n^oo a^_-^(2u;o,c) 

exists, which explains the second numerical observation of section [lT| 



3.5 Theoretical explanation to the third numerical observation 



In sections |3.3| and 3.4 we focussed the discussion on the asymptotic behavior for families satis- 
fying hypothesis H2'. The map considered in the third numerical observation of section [LT] is 
an example of a map not satisfying H2' for which equations (33) and (35) do not hold. 



Let {c(a, e)}(^ be a two parametric family of maps satisfying hypothesis HI, H2 and 
H3. Let a^(w,c) denote slope of the reducibility loss bifurcation associated to the 2" periodic 
invariant curve of the family. Finally consider ujq a Diophantine rotation number for the family. 
Let a* be the parameter value for which {c(a, 0)}(Q^o)eA intersects W'^{(^,1Z). 

The main difference between the example considered in the third numerical observation of section 



1.1 and the previous two is that 

dec{a*,Q) = vo^i + fo,2 with vq^i € Bi, i 



1,2, 



where the spaces Bi are given by (10). More concretely, for the numerical example cited above 
we have considered 

vo,i = fi{x) cos{9), vo,2 = nh{x) cos(26'). (36) 
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As this family depends on rj, we denote by Crj this concrete family. The parameter rj is considered 
in addition to the parameters a and £ of the family. In other words, for each > 0, is a two 
parametric family of maps. 



Then one has that equation (|28|) still holds with Vn = Vn,i + fn,2, with Vn,i and Vn,2 the sequence 

In this case there is no universal behavior because 
where the renormalization operator is not 



given by (34) and fo,i and fo,2 given by (36) 

Vn lives in a bigger invariant space 



the sequence 

contractive. On the other hand, the numerical computations in [21j suggest that the sequence 
a'n{u}o, Cr^) / a'^_i{ujQ, Cfj) (for 77 > 0) is not asymptotically equivalent to a^(t<;o, cq )/«'„_ ;^(ti;0) cq), 
but both sequences are r/-close to each other. This can be explained as a consequence conjecture 
D where we conjecture uniform growth (in norm) of the sequences Vn,i and Vn.2- For more 



details see theorem 3.20 in [23j. 



4 Applicability to the Forced Logistic Map 

The theory exposed in sections [2] and [3] have been built as a response to the observations 
done in the study of the Forced Logistic Map (see |1H I21j). In this section we discuss the 
applicability of the quasi-periodic renormalization theory to the Forced Logistic Map. In the 
cited papers we considered two different version of the FLM, which correspond to the map 
([2]) with either ha,£{0,x) = ax{l — x) cos(27r0) or ha,£{9,x) = cos(27r6'), where the parameters 
(a, e) G [0, 4] x [0, 1]. Notice that these two forms of the FLM do not satisfy the requirements of 
the theory developed in the previous sections because the family of maps does not belong to B. 

This problem can be easily solved as follows. For a > 2 we can consider the affine change of 
variables given hy y = ax + b, with a = and h = — ^j^- If we apply this change of variables 
to the family ([2]) when ha,£{0,x) = ax{l — x)cos(27r^) we obtain the following family, 

e = e + io, 

fa a{a- 2) A , ,0 2 V (37) 

y = a (l + ecos(27r6')) ' ^ ^ 



a — z 4 / a 



If we apply the same change of variables when ha^ei^^j ^) — cos(27r0) we obtain this other family 

e = e + io, 

«(«-2) 2 4e (38) 

y = 1 y + -cos27r6'. 

4 a — 2 

Although the change of variables considered depends on the parameter a, the parameter space of 
the maps ( [37[ ) and ( 38 ) is the same as the parameter space of the map ^ (for the corresponding 



value of h). Then any conclusion drawn on the parameter space of the map (37) (respectively 
(38)) extends automatically to the parameter space of ([2]). 



With this new set up, we have that both families of maps belong to B for a G (2, 4) and e 
small enough. One should check that the FLM satisfies hypotheses HI and H2'. To check HI 
one should check that the one dimensional Logistic Map intersects transversally W^{^,Tl) the 
stable manifold of the renormalization operator. This is an implicit assumption when one uses 
the renormalization operator to explain the universality observed for the Logistic Map. The 
only proof (to our knowledge) of this fact is the one given by Lyubich (theorem 4.11 of [17J). 
This proof is done in the space of quadratic- like germs, which is a smaller space than the one 



considered here. Hypothesis H2' is trivial to check for the maps (37) and (38)). 



23 



n 




ea 


e-r 


1 


-5 


832915e+00 


1 


986668e- 


15 


3 


405962e- 


16 


2 


-8 


494260e+00 


1 


155106e- 


13 


1 


359866e- 


14 


3 


-1 


635128e+01 


1 


265556e- 


14 


7 


739800e- 


16 


4 


-1 


125246e+01 


3 


995769e- 


14 


3 


551018e- 
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224333e+01 
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009061e- 


13 
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-1 


807969e+01 
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971921e- 


13 
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-3 


473523e+01 
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087996e- 
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328470e- 


12 
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-2 


958331e+01 
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204433e- 
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7 


451609e- 


12 


9 


-4 


156946e+01 


1 


293339e- 


09 


3 


111273e- 


11 


10 


-7 


896550e+01 


3 


847537e- 


08 


4 


872428e- 


10 


11 


-7 


450073e+01 


7 


779131e- 


08 


1 


044168e- 


09 



Table 2: Values of Oi'^{uj) for the FLM map (37) and uj = The values £a and correspond 



to the discrepancy with the estimates of |21j in absolute and relative terms. 



On the other hand, note that theorem 3.8 in [22j not only gives the existence of reducibility loss 
bifurcations, but it also gives an explicit expression of its slopes in term of the renormalization 



operator (the ones given by formulas ( 18 ) and ( 19 )). Actually we have given even a more explicit 



formula in corollary 3.13 in |22j . We used these formulas to compute the reducibility directions 



a'n{ijj) of the Forced Logistic Map (37) 



The initial values q„ have been computed numerically, by means of a Newton method applied to 



their invariance equation. Differentiating on formula ( |37[ ) (respectively (38)) it is easy to write 
the values of Uq"^ and Uq"^ in terms of an- Then using the discretization of the operator done 

p(n) 



a' 
we can compute numerically the iterates f^\u^k^ and for /c = 1, n — 1. 
Once we have these functions, we can evaluate them to compute the values of a^(w) given by 
formula (48) of 



We have used this procedure to compute the values of a^(f^) for the map (37) (consequently 
also for the map (2) with ha,e{G,x) = ax{l — x)cos{2'k6)) and = . The results are 
shown in table [2| The values 0^(0;) have been also computed in |21j via a completely different 
procedure, based on a continuation method with extended precision. More concretely these 
values are displayed in table 2 of the cited paper. The third and fourth columns of table [2] 
display the discrepancies between both computation, in absolute (ea) and relative (e,.) terms. 



This experiment has been repeated for other values of uj and also for the map ( 38 ) . In all cases 
we obtained that the slopes computed with both methods are the same up to similar accuracies 
to the ones displayed in table [2] 



This supports the correctness of both computations and at the same time the conjecture A 
which has been assumed to be true to derive the formula used for this estimation. 



A Numerical computation of the spectrum of £^ 

In this appendix we present the numerical method that we used to discretize C^j and to study 
its spectrum numerically. Our method is a slight modification of the one introduced by Lanford 
in [15] (see also [16j). 
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Let D{zq,p) be the complex disc centered on zq with radius p. Consider Tn-l(D{zQ, p)) the 
space of real analytic functions such that they are holomorphic on D(zQ,p) and continuous on 
it closure. Given a function / G TZ'H{D{zo, p)), we can consider the following modified Taylor 
expansion of / around zq, 

fc=0 \ / 

The truncation of a Taylor series at order N induces a projection defined as 

/ ^ (/o,/l, • • • ,/Af), 

On the other hand we have its pseudo-inverse by the left 

i(^) : M^+i ^ 71^ (Dp) 

iV / _ \i 

(/o,/i,...,/^) ^ E/^^J ' 

in the sense that i(7v) ° P{n) is the identity on M^"*""^. Note also that both maps are linear. 

Given a map T : TZ'H{D{zo, p)) — >• TZ'H{D{zq, p)), we can approximate it by its discretization 
T(^) : M^+^ — M^+-'^ defined as T^^^ := P(Ar) oToi^^-j. More concretely if T is a linear bounded 
operator, we can compute the eigenvalues of the operator T^^^ in order to study the spectrum 
of T. In general the eigenvalues of T^^-* might have nothing to do with the spectrum of T. 
For example an infinite-dimensional operator does not need to have eigenvalues, but a finite- 
dimensional one will always have the same number of eigenvalues (counted with multiplicity) as 
the dimension of the space. For this reason will do some numerical test on the results obtained 
with this discretization. 

At this point consider the map : 7^'H(W)e7^'W(W) 7^'H(W)e7^'W(W) defined by equation 
(11). If we set W = D{zo, p) we can use the method described above in each component of 
nniW) e nniW) to dlscretize and approximate it by a map cl^'' : M^CA'+i) _^ m2(7V+i)_ 
Concretely in our computation we have taken zq = g and p = | . In figure [2] we include graphical 

" " can be found more details about 



2.1 



evidence that the set W = D (g, |) satisfies HO, in section 
this. 

In table [3] we have the first 24 eigenvalues of for N = 100 and uj = . The eigenvalues 
have been sorted by their modulus, from bigger to smaller. Note that the eigenvalues of the 



discretized operator also satisfy the properties given in section 2.1 To justify the validity of 
these eigenvalues we have done the following numerical tests. 

Consider that we have a real eigenvalue of multiplicity two, or a pair of complex eigenvalue^ 
which is persistent for different values of N (the order of the discretization). The first test 
done to the eigenvalues is to check if the distance between the associated eigenvectors decreases 



when is increased. In the left part of figure 10 we have the distance between the eigenvectors 



associated to the same eigenvalue of the operators ci^^ and cj^^^'^ as a graph of N, with 
varying from 40 to 100. We have plotted this distance for the first twenty-four eigenvalues. 
To compute the distance between eigenvectors we have estimated the supremum norm of the 
difference between the real function represented by each of the vectors, in other words we have 
computed \\i'^^\v^^^) - z(^i°)(7;("°))||oo in the interval {zq- p,zo + p) = WnM. 



^Note that there are no simple eigenvalues due to corollary 2.18 in [22] 
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0000000? 
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+0 
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+0 


0869398f 
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+0 
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Table 3: The first twenty four eigenvalues of CIj , for ^ = — 2 • ^'^^ computation N has 
been taken equal to 100. 
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0.0001 *\ 
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Figure 10: Estimation of the errors and the radii of convergence of the first twenty-four eigen- 
values of for oj = with respect to the order of the discretization. See the text for more 
details. 



Note that, since the distance goes to zero this indicates that the eigenvectors, namely v^^\ 
converge to a limit v*. One should expect these eigenvalues to be in the spectrum of jC^j, but 
nothing ensures that v* belongs to the domain of jOi^. We have done a second test on the 
reliability of the approximated eigenvectors, where we check this condition. 

Let us remark that with the numerical computations done so far, we have only checked that the 
eigenvectors as elements of 7^'H(D(zo, p)) ©7^?^(D(zo, p)) to converge in the segment {zq — p,zq + 
p) C M instead of checking the convergence in the whole set D{zq, p). Let us give evidences that 
approximate eigenvectors obtained with our computations have a domain of analicity containing 
O{zo,p). 

Consider that we have a function / holomorphic in a domain of the complex plane containing 



D(zo, p). If the we consider the expansion of / given by equation (39 ), we have that r the radius 
of convergence of the series around zq is given as 

^ ^ P 

lim sup„_5.(^ 
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Figure 11: Estimation of the distance between eigenvectors with the same eigenvalue (left) and 
estimation of the radios of convergence (right) of the first twenty-four eigenvalues of Ci^ with 
respect to to. See the text for more details. 



With the discretization considered here we have an approximation of the terms fn, hence these 
can be used to compute a numerical estimation of r. 



Consider v an eigenvector of the operator Cui- We have that v = {vi,V2) £ TZ7i{W) (B'R'H{W). 
Given = {v['^\v2^^) a numerical approximation of the eigenvector, we can use the procedure 
described above to estimate the radius of convergence of each vi and V2. We have done this 
for the eigenvectors associated to each of the first twenty-four eigenvalues of C^j with to = 
(keeping only the smaller of the two radius obtained) . The results are displayed on the right part 
of figure \T0\ where the estimated radius has been plotted with respect to N, the order of the 
discretization. Note that the estimations give a radius bigger than p 
the eigenvectors are analytic in 



|, which indicates that 
zo,p), and continuous on its closure, for zq = I and p = §• 



Up to this point, we have considered w fixed to but the same computations can be done 

to study the spectrum of C^i with respect to the parameter u. In figure [3] of section 2J_ we 
have plotted the first twenty- four eigenvalues of the map with respect to uj. The set T has been 
discretized in a equispaced grid of 1280 points. Recall that the operator depends analytically 
on CO (proposition 2.20 in [22]), therefore the spectrum also does (as long as the eigenvalues do 
not collide, see theorems III-6.17 and VII- 1.7 in [13j). The numerical results agree with this 
analytic dependence. 

For this computation we have also made the same test as before to the eigenvalues. The results 



of these tests are shown in figure 11 To estimate the convergence of the eigenvectors we have 
compared the eigenspaces of the eigenvalues of C^^^ with the eigenspaces associated to cl^^^^ 
for each value of u in the cited grid of points on T. The estimation of the radius of convergence 
has been also done with respect to uj for = 90. We have plotted the estimated error and 
convergence radius for the first twenty-four eigenvalues in the same figure. Both result indicate 
that the eigenvalues obtained are reliable. 
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